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The magnetic hyperfine field and electric-field gradient at isolated lanthanide impurities in an Fe 
host lattice are calculated from first principles, allowing for the first time a qualitative and quantita- 
tive understanding of an experimental data set collected over the past 40 years. It is demonstrated 
that the common Local Density Approximation leads to quantitatively and qualitatively wrong re- 
sults, while the LDA-fU method performs much better. In order to avoid pitfalls inherent to the 
LDA-I-U method, a careful strategy had to be used, which will be described in detail. The lanthanide 
4f spin moment is found to couple antiferromagnetically to the magnetization of the Fe lattice, in 
agreement with the model of Campbell and Brooks. There is strong evidence for a delocalization- 
localization transition that is shifted from Ce to at least Pr and maybe further up to Sm. This shift 
is interpreted in terms of the effective pressure felt by lanthanides in Fe. Implications for resolving 
ambiguities in the determination of delocalization in pure lanthanide metals under pressure are dis- 
cussed. For the localized lanthanides, Yb is shown to be divalent in this host lattice, while all others 
are trivalent (including Eu, the case of Tm is undecided). The completely filled and well-bound 5p 
shell of the lanthanides is shown to have a major and unexpected infiuence on the dipolar hyperfine 
field and on the electric-field gradient, a feature that can be explained by their dependence. 
An extrapolation to actinides suggests that the same is true for the actinide 6p shell. The case of 
free lanthanide atoms is discussed as well. 

PACS numbers: 71.20.Eh, 75.20.Hr, 75.25.+Z, 76.80.-|-y, 76.60.-k 



I. INTRODUCTION 



A prototype problem in the field of nuclear condensed 
matter physics is to determine and understand the mag- 
netic hyperfine field (HFF) at any of the elements of the 
periodic table, incorporated as a substitutional impurity 
in a simple ferromagnetic host such as bcc Fe. Under- 
standing hyperfine fields forms a critical test for our un- 
derstanding of condensed matter. Moreover, they pro- 
vide a convenient tool for nuclear physicists to determine 
nuclear magnetic moments: two features that explain the 
decades of experimentali and theoretical^i^ efforts that 
have been devoted to this problem. Today, the hyperfine 
fields of all elements as substitutional impurities in bcc 
Fe are well-understood up to about Z=55^i^*SiL&. For the 
heavier 5d impurities, sizeable deviations between theory 
and experiment remain^. The hyperfine fields of very 
light impurities at interstitial sites in Fe have been calcu- 
lated as welUS*ii. Lanthanide impurities in Fe are much 
less understood, both experimentally and theoretically. 
As far as experiment is concerned, it is hard to obtain 
reliable values for the lanthanide hyperfine fields (and 
also for their electric-field gradients, see below). This 
problem is illustrated by the rather desperate conclusion 
of L. Niesen in a still useful reviewi^ back in 1976: 
if we cannot perform experiments that yield unambiguous 
results, we should better do no experiments at alir On the 
theoretical side, no ab initio studies have been performed 
yet for lanthanide impurities in Fe (an approach using a 
model Hamiltonian is developed in Ref. Il3() . The reason 
for this lies in a known failure of the widely used Lo- 
cal Density Approximation (LDA) within Density Func- 
tional Theory (DFT): LDA is not suitable to describe 



strong electron correlationsi^iiSiiSiii. As a result, the 
strongly correlated and mainly localized 4f states in lan- 
thanides are rendered itinerant by LDA. It can therefore 
be anticipated that for a lanthanide impurity in bcc Fe, 
LDA is incapable to describe correctly the interaction 
between these localized and strongly correlated 4f elec- 
trons and the itinerant 3d states of the host material. 
This could be overcome by treating the 4f states as core 
states, an approach that has been used in the past with 
some success for lanthanide9i^*i^*ii. In this way, how- 
ever, one forces the f-electrons to behave exactly as in 
free atoms, which is not entirely correct. An efficient 
and popular way to improve on the LDA failure with- 
out resorting to fully atomic 4f behavior, is to use the 
LDA-f U methodi2*i2iSL2ii2a. In LDA-hU, the correlation 
absent in LDA is reintroduced by an on-site Coulomb re- 
pulsion parameter U, to which an a priori value has to 
be assigned. The LDA-I-U method has been used in the 
recent past with considerable success (recent examples 
are Refs. I2lll2l2l2l2^ and many others), but it is not 
yet clear where the boundaries of its range of applicabil- 
ity are. In this work, we will examine how well LDA+U 
performs on a delicate quantity as the HFF. 

An extra feature for lanthanide impurities in Fe that 
is absent for lighter impurities, is the presence of a large 
electric-field gradient (EFG) at the lanthanide nucleus. 
At a site with cubic point symmetry - such as a sub- 
stitutional site in bcc Fe - the EFG tensor must be 
necessarily zero. Due to the strong spin-orbit coupling 
for these heavy impurities, however, the crystalline cu- 
bic point symmetry at the lanthanide site is lowered to a 
tetragonal one. This allows the existence of a large EFG 
at the nucleus of the impurity. The same happens for 
5d impurities in Fe2L22iSS, but there the EFG is 2 orders 
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of magnitude smaller. We have calculated and analyzed 
this EFG for lanthanides, and compare it with the sparse 
experimental data. 

The goals of this work can be summarized as follows. 
On the physical side, we want to obtain better quantita- 
tive and qualitative insight in magnetic hyperfine fields 
and electric-field gradients of lanthanide impurities in Fe. 
This should allow us to asses better the reliability of the 
existing experiments, and to derive the underlying phys- 
ical mechanism. On the technical side, we want to ex- 
amine whether the range of applicability of the LDA-f U 
method can be extended to problems as delicate and sen- 
sitive as magnetic and electric hyperfine interactions of 
heavy impurities in a transition metal host. It will be 
shown that in the course of this analysis unexpected new 
results and questions show up, such as the influence of the 
lanthanide 5p electrons fSec. lIII Cl and llVB() and the po- 
sition of the delocalization-localization transition in this 
system (Sec. ImTl and rVll) . 



II. COMPUTATIONAL DETAILS 

All our calculations were performed within Density 
Functional Theorji2Si2ii2S, using the Augmented Plane 
Waves -I- local orbitals (APW-f lo) method^SiSi^ as im- 
plemented in the WIEN2k package^^ to solve the scalar- 
relativistic Kohn-Sham equations. In the APW-flo 
method, the wave functions are expanded in spherical 
harmonics inside nonoverlapping atomic spheres of ra- 
dius i?MT, and in plane waves in the remaining space of 
the unit cell (=the interstitial region). For the Fe atoms 
a Rmt value of 2.20 a.u. was chosen, while for the lan- 
thanide impurity we used Rmt =2.45 a.u. The maximum 
£ for the expansion of the wave function in spherical har- 
monics inside the spheres was taken to be £max = 10. 
The plane wave expansion of the wave function in the 
interstitial region was made up to i^max = 7.5/ R^^^ = 
3.41 a.u.~^, and the charge density was Fourier expanded 
up to Gmax = IGVRy- 

The lattice constant of Fe was fixed at the experimen- 
tal value of 2.87 A. In order to reproduce the situation 
of an isolated impurity in bulk Fe, we used the super- 
cell approach with a 2 x 2 x 2 supercell where one iron 
atom was replaced by a lanthanide atom. The neigh- 
boring Fe atoms will be displaced by the presence of this 
impurity, as was documented before for lighter impurities 
in Fa^i^. We took this effect into account in an average 
way by relaxing the nearest neighbors for Eu as an impu- 
rity (which is in the middle of the lanthanide series) , and 
kept the same relaxation fixed for all other lanthanides. 
The Eu-Fe distance was 2.60 A which is an increase of 
0.11 A with respect to the Fe-Fe distance and which is 
almost identical to the distance between 5p impurities 
and their Fe neighbors^. It was tested for another lan- 
thanide (Er) that there was only a marginal difference of 
less than 1 T between the HFF obtained with the Eu-Fe 
distance and the correct Er-Fc distance (2.58 A). A test 



for an extended supercell of 32 atoms was also performed. 
We relaxed the first four nearest neighbors. The Eu-Fe 
distance hardly changed (2.63 A) and the Fermi contri- 
bution to the hyperfine field changed with 5 T. For the 
sampling of the Brillouin zone (BZ) of the 2x2x2 su- 
percell we took 75 special k-points in the irreducible part 
of the BZ, which corresponds to a 10 x 10 x 10 mesh. 

As exchange-correlation functional, the Local Density 
Approximation (LDA)^ was used. Spin-orbit (SO) cou- 
pling was taken into account in all the calculations by 
a second-variational step scheme^, using a cut-off en- 
ergy £;fO = 3.0 Ry. Relativistic Local Orbitals (RLO) 
for the lanthanide 5p states were added to the basis set, 
because it is known that for actinides this allows to re- 
duce the basis set size needed for the second variational 
stepS (=lower i?f^). Limitations in the implementa- 
tion prevent to obtain correct EFG's and dipolar HFF's 
when RLO's are used. Therefore, whenever such infor- 
mation was needed, the calculations were repeated with- 
out RLO's. This never had a large influence on the ob- 
tained values, however. For the LDA-I-U method, the 
'Around the Mean Field' (AMF) scheme of Czyzyk and 
SawatzkjiiS was used. The choice of the U and J param- 
eters is discussed in detail in Sec. IIII CI 

Free atoms were simulated by a supercell containing 
only one lanthanide and vacuum otherwise, leading to 
a separation of 9.4 A between two 'neighboring' lan- 
thanides. For free ion calculations, this cell was charged. 
All the other parameters were chosen exactly the same 
as in the calculations for lanthanides in Fe. 



III. MAGNETIC HYPERFINE FIELDS 

A. Experimental data set 

Let us first have a look at the experimental data set for 
the HFF (Fig^Ja). Only in 4 cases the magnitude of the 
HFF and its sign are known with high reliability (the sign 
of the HFF indicates whether the field is parallel (+) or 
antiparallel (-) with respect to the magnetization of the 
Fe host lattice). These cases are La (-47(1) T)^, Ce 
(-41(2) T)ia and Lu (-73.12(36) T)ii for which Nuclear 
Magnetic Resonance on Oriented Nuclei (NMR/ON) has 
been performed, and Yb (-125(8) T)i^ on which Time- 
Dependent Perturbed Angular Correlation spectroscopy 
(TDPAC) has been applied. The latter technique has 
also been used for Gd^, albeit on a recoil-implanted 
sample which is not necessarily clean. The value of - 
26(8) T obtained in this way agrees well with an in prin- 
ciple reliable Mossbauer measurement of -37 T, which is 
unfortunately not very well documented^. Three time- 
integrated Perturbed Angular Correlation (IPAC) mea- 
surements are available for Gd as well - IPAC is a method 
that is rather unreliable, and can merely be used to 
determine the sign and an order of magnitude. They 
yield -20(5) Ti^, -18(9) Ti& and -7 Til. A HFF of - 
30(10) T can therefore be assigned to Gd in Fe in a reli- 
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FIG. 1: a) Experimental data set for the magnetic hyperfine 
fields of lanthanides in Fe. If the sign of the HFF is not mea- 
sured, the data point is plotted both at positive and negative 
values. Distinction is made between highly reliable data for 
which the sign is measured (diamond), highly reliable data 
without sign measurement (square), less reliable data with- 
out sign (triangle) and data that are rather unreliable for the 
magnitude of the HFF but reliable for the sign (circle). For 
references and values, see text. The line connects the most 
likely values for all lanthanides. If multiple measurements 
with the same reliability were available, only one of them is 
given. More data can be found in the compilation of Rao-, 
b) Comparison between experiment and several types of cal- 
culations for the magnetic hyperfine field of lanthanides in Fe. 
Diamonds: most probable experimental data points (this is 
the full line from Fig.0-a). Dotted-dashed line: LDA results 
for the antiferromagnetic orientation (Campbell-Brooks orien- 
tation). Dotted line: LDA results for the ferromagnetic orien- 
tation when this one has the lowest energy. Full line: LDA-I-U 
value for trivalent lanthanides. Dashed line: LDA-fU value 
for divalent lanthanides. 



able way. In 4 other cases the magnitude of the HFF but 
not its sign has been measured with an accurate method 
as Mossbaucr Spectroscopy (MS): Eu (148.2(9) T)i24a, 
Dy (610(7) T)a, Er (768(13) T)^ and Tm (671 T>i^. 
For Pr— and Tb^, the magnitude but not the sign has 
been measured with Low Temperature Nuclear Orienta- 
tion (LTNO). This non- resonant technique provides data 
that are less accurate than the previous ones, although 
they still are reasonably reliable. Finally, in the case 
of Co^, Nd^, SmS, Eu^, Gd^, Dj4&, Eri^ and TuM. 
IPAC experiments have been reported, from which only 
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FIG. 2: Schematic summary of the model of Campbell and 
Brooks. The localized 4f spin moment does not interact di- 
rectly with the moments of neighboring atoms. But it does in- 
teract with the 5d moment on the same atom by ferromagnetic 
intra-atomic exchange. The more delocalized 5d spin moment 
interacts by antiferromagnetic interatomic exchange with the 
Fe 3d moment. As a result, the net interaction between Fe- 
3d and lanthanide-4f is antiferromagnetic. A consequence of 
the model of Campbell and Brooks is the positive (negative) 
HFF in the first (second) half of the lanthanide-series in Fe 
(Fig. 0a). Due to Hund's third rule, the lanthanide 4f orbital 
moment is antiparallel (parallel) to the 4f spin moment in the 
first (second) half of the series. Because the orbital contribu- 
tion to the HFF (which is parallel to the orbital moment) is 
dominant for lanthanides, the total HFF is parallel (antipar- 
allel) to the Fe moment - and hence called positive (negative) 
- in the first (second) half of the series. 



the sign information can be reasonably trusted (see e.g. 
the agreement with other experiments in Fig. ^a and 
Ref. [l|) . Due to the the latter sign information, the HFF 
of the light lanthanides is guessed to be positive, while for 
the heavy lanthanides it is negative. The line in Fig. 
a summarizes the most likely interpretation of this data 
set. 

Fig^a can be understood in terms of Hund's rules and 
the model of Campbell and Brooks. Based on heuris- 
tic arguments (Ref. Is^) and first principles calculations 
(Ref. l59|) ■ Campbell and Brooks showed that the inter- 
atomic exchange interaction between a transition metal 
3d spin moment and a lanthanide 5d spin moment is 
antiferromagnetic (Fig. [5J). The lanthanide 4f moment 
is localized at the lanthanide site and cannot directly 
interact with its transition metal neighbors, but it has 
a ferromagnetic intra-atomic exchange interaction with 
the lanthanide 5d moment. The result is a net antifer- 
romagnetic coupling between the lanthanide 4f moment 
and the transition metal 3d moment (Fig. Accord- 
ing to Hund's third rule, the lanthanide orbital moment 
is antiparallel to the lanthanide spin moment for the 7 
lightest lanthanides, and parallel to it for the 7 heaviest 
lanthanides. The dominant contribution to the HFF is 
the orbital HFF (see Sec lHIBl and Fig.EJ, which IS par- 
allel to the orbital moment. Therefore, one expects the 
total HFF to be parallel to the Fe magnetization (and 
hence positive) for the light lanthanides, and antiparallel 
(negative) for the heavy ones, as is seen indeed in Fig.^ 
a. 
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FIG. 3: White circles indicate the ferromagnetic solution 
(FM), black triangles the antiferromagnetic solution (AF). 
Symbols are LDA results, lines are LDA-f U results (obtained 
from Hund's rules occupations, see text). Mind the scale, 
which is 10 times larger in a) compared to b) and c). a) Or- 
bital contribution to the HFF at the lanthanide site due to 4f 
electrons, b) Dipolar contribution due to 4f only (fully line), 
and due to 4f and 5p (dashed line, see text), c) Fermi contri- 
bution. The total HFF is the sum of these 3, and is almost 
undistinguishable from a). 



B. LDA calculations 

As a first step, we calculate the magnetic HFF with 
the common LDA. This will provide us with a data set 
to which we can later compare the possible improvement 
by LDA-I-U, and it allows to introduce some pccuUaritics 
that will play a role in all later calculations as well. As is 
usual with this type of methods, our calculations involve 
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FIG. 4: Energy difference between the situations with the 
lanthanide 4f spin moment ferromagnetically aligned with the 
Fe 3d moment and with the 4f moment antiferromagnetically 
aligned (Efm - Eaf). If this energy difference is positive, 
the antiferromagnetic situation has the lowest energy. En- 
ergy difference for LDA are compared with energy differences 
for LDA-I-U with progressively larger U. The LDA result is 
equivalent to U=0.0 Ry. There are some gaps in the picture, 
because not for every situation a converged solution could be 
found. For Pm the calculated energy difference is much larger 
(0.42 Ry for U=0.6 Ry) and is out of scale. 



an iterative procedure ('sclf-consistcnt field' procedure) 
that yields in the end a possible state of the calculated 
system, which is not necessarily the desired ground state: 
in the space of possible solutions, this self-consistent field 
procedure finds a local minimum, but not necessarily the 
global minimum. The local minimum that is obtained, 
depends to some degree on the starting configuration 
that was initially chosen. This behavior is prominently 
present for lanthanides in Fe. If the spin moment of the 
lanthanide initially is put parallel to the Fe spin moment, 
then this parallelism is maintained throughout the iter- 
ative procedure (except for La, Ce and Pr, where the 
moment always spontaneously turns to an antiparallel 
orientation). We call this from now on the ferromagnetic 
solution. With an initially antiparallel configuration, an 
antiparallel (or antiferromagnetic) solution is obtained. 
If the lanthanide was given initially no spin moment, then 
a solution with a spin moment that is much reduced com- 
pared to the two preceding solutions was found. In order 
to decide which of those is the ground state, one has to 
look at the total energy of each solution. The total en- 
ergy of the case with reduced moment was much higher 
than the others, and we will not consider it further. The 
energy differences between the other two solutions are 
given in Fig. 0] For all lanthanides up to Tb, the an- 
tiferromagnetic solution has the lower energy. Starting 
with Dy, the ferromagnetic solution becomes the ground 
state. 

In Fig. 121 the different contributions to the magnetic 
HFF arc given for both types of solutions. A HFF is 
a magnetic field at the position of the nucleus, and it 
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is built mainly from 3 contributions: the spin dipolar 
field, the Fermi contact field and the orbital field. The 
spin dipolar field is generated by the spin moments of 
the electrons surrounding the nucleus. For cubic point 
symmetry, this contribution vanishes. The Fermi contact 
fieldSS is of dipolar nature as well, but is due to the pene- 
tration of s-electrons into the nucleus. It does not vanish 
for cubic symmetry, and it is the dominant (and almost 
only) contribution for impurities up to Z=55 in Fe. The 
orbital field stems from the electric charge of the electron 
that orbits the nucleus, and it vanishes for cubic symme- 
try. As Fig. 121 shows, the dipolar field does not exceed 
a few Tesla, while the Fermi contact field lies between 
-20 T and -40 T. The orbital field is the dominant contri- 
bution, and can reach almost ±600 T. At first sight, one 
would expect a zero orbital and dipolar field for a substi- 
tutional impurity in Fe, as the point symmetry is cubic. 
The reason why this is not the case for lanthanides, is 
purely due to the presence of spin-orbit coupling, which 
breaks the crystalline cubic point group symmetry. The 
oscillatory behavior of the orbital moment reflects Hund's 
third rule: the orbital moment is antiparallcl (parallel) 
to the spin moment in the first (second) half of the lan- 
thanide series. Because the orbital field is parallel to the 
orbital moment, it will for the antiferromagnetic solution 
be positive in the first half of the series and negative 
in the second half (and vice versa for the ferromagnetic 
solution). 

According to the LDA total energies, we have to accept 
the antiferromagnetic solution as the ground state up 
to Tb, and the ferromagnetic solution starting from Dy. 
This leads to positive HFF's for almost all lanthanides 
(Fig. Q-b), which is contradiction with the current in- 
terpretation of the experimental data set and with the 
model of Campbell and Brooks. We will demonstrate in 
Sec. 1111 CI that this is a new, clear example of a failure 
of LDA, to be added to the list of notorious shortcom- 
ings as the general overbinding behavior, the prediction 
of the wrong crystal structure for Fe, and the prediction 
of metallicity for some strongly correlated insulators as 
NiO and La2Cu04. Even if we would select the antiferro- 
magnetic solution throughout (as the Campbell-Brooks 
model suggests), then still the quantitative agreement 
with the experimental HFF's is rather poor (Fig.^b). 



C. LDA+U calculations 

Using and interpreting LDA-f-U calculations brings 
some complications that are absent for LDA. First, 
LDA+U schemes are not fully ab initio: they involve 
an on-site Coulomb repulsion parameter U and an on- 
site exchange interaction constant J that have to be cho- 
sen a priori for every atom with strong correlations. In 
our case we have to choose one U and one J for the f- 
states of the lanthanide impurity. In line with the strat- 
egy adopted for the relaxation (Sec. we strive for 
reasonable overall agreement and do not focus on agree- 



ment for individual cases too much. Therefore we take 
the same U and J for all lanthanides. Looking at other 
calculationsiSiSi and experimentsSSiS, JJ = 0.6 Ry is a 
reasonable choice. The value of J is usually an order of 
magnitude smaller, and it does not affect the results as 
much as U does. Therefore wc take J as 10% of U. An- 
other - second - complication with LDA-f-U is that due 
to the introduction of tuneable variables U and J, the 
total energy looses its exact mathematical meaning as a 
variational quantity. If LDA-I-U is used, a lower energy 
for one solution does therefor not necessarily mean that 
this solution is to be preferred over another one with a 
higher energy. This does not mean that these total en- 
ergies are meaningless: for instance, LDA-I-U has been 
used successfully to determine the equilibrium volume of 
materials by energy minimization^. The current gen- 
eral opinion is that if the electronic structure does not 
qualitatively change between two slightly different situ- 
ations, then the total energies obtained by LDA+U for 
these situations can probably be compared. The crite- 
rion to decide whether there is a qualitative difference 
between two electronic structures or not, is provided by 
the 4f density matrix: in qualitatively similar cases, the 
occupation of the 14 m-orbitals (= the diagonal elements 
of the two spin-polarized 7x7 density matrices) should 
be more or less identical. Finally - and third - there are 
much more local minima in the space of solutions when 
LDA+U is used, and a calculation gets easily trapped in 
one of them. Because the density matrices for two dif- 
ferent solutions are necessarily different from each other, 
the total energy cannot be used to determine which so- 
lution is the ground state. This problem is illustrated in 
Tab. n where the diagonal elements of the f-up density 
matrix is given for various antiferromagnetic solutions for 
Tm in Fe. The 5 electrons can be distributed in different 
ways over the 7 orbitals, and always a converged solution 
can be obtained. The HFF field can be very different for 
all cases. In the second column, the total energy of these 
5 solutions is given, relative to the case with the lowest 
energy (case 2). Case 1 - which we will later identify as 
the most probable ground state - is only third in rank if 
the total energy is considered: an illustration of the lack 
of meaning of the LDA+U total energy for cases with 
different density matrices. This makes it impossible to 
proceed in a straightforward way as was done for LDA. 

Faced with this problem, we first take one step back 
and examine the energy difference between the ferromag- 
netic and the antiferromagnetic case for the same type 
of solution (i.e. for a ferro- and and antiferromagnetic 
solution that have a similar f-electron density matrix). 
These energy differences can be expected to be meaning- 
ful (see above), a hope which is supported by the fact 
that we find them not to depend on the type of solution 
for which we make the comparison. The result is shown 
in Fig.H for 3 different values of J7: 0.2, 0.4 and 0.6 Ry 
(some gaps arc present, because not for every lanthanide 
a solution could be found for every U). Clearly, the use 
of a non-zero U makes the antiferromagnetic case more 
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TABLE I: Diagonal elements of the 7x7 4f-up density matrix for Tm in Fe (antiferromagnetic case), together with the orbital 
(4f), dipolar (4f+5p) and Fermi contributions to the total HFF (Tesla). These diagonal elements give the occupation of each 
m-orbital (between and 1). The LDA result is compared with several LDA+U calculations, all with U=0.6 Ry. The LDA+U 
calculations differ only in the initial distribution of the f-electrons over the different orbitals. In the second column, the total 
energy (in mRy/atom) of the LDA+U calculations is given, relative to case 2, which has the lowest energy (see text for 
discussion) . 
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stable, for all lanthanides. This brings the sign of the 
HFF in agreement with experiment and with the model 
of Campbell and Brooks. We conclude that LDA+U de- 
scribes the effective d-f exchange interaction much better 
than LDA does, and that LDA is qualitatively wrong in 
this respect. 

Next, we try to find a way to obtain ground state val- 
ues for the HFF in the antiferromagnetic case. To this 
end, we turn the annoying freedom of having several ways 
to occupy the m-orbitals into an advantage, by determin- 
ing the individual contribution of each m-orbital to e.g. 
the orbital field. This can be done by first calculating 
the orbital HFF for several different antiferromagnetic 
solutions, as is given as an example for Tm in Tab. 
Then a system of linear equations is set up, with as 7 
variables Xm the orbital HFF of each of the 7 m-orbitals. 
The occupation of each of these orbitals (or the diagonal 
elements of the density matrix from Tab.|3) are the coef- 
ficients. The occupation found in the calculation should 
give the calculated orbital field, as is illustrated here for 
'case 1' in Tab. HI 

1.00a;_3 + 0.99a;_2 + 0.99x_i + 0.99x0+ 

0.99a;i + 0.01x2 + 0.01x3 = -718 

This system of equations can be supplemented by other 
equations expressing some general truths (the orbital 
field with all 7 m-orbitals filled is zero, the contribution 
by +m is opposite to the one by -m), such that the system 
becomes overdetermined. Each subset of 7 independent 
equations should give the same x^, which indeed they 
do. 

This was done for 10 elements from the lanthanides 
series (Ce, Nd, Pm, Sm, Tb, Dy, Ho, Er, Tm and Yb), not 
only for the orbital HFF but also for the orbital moment 
and the dipolar HFF. The dominant contribution to the 
orbital HFF is due to the 4f electrons, as one could expect 
(Fig. |51-a,b) . Orbitals with opposite m quantum number 
yield opposite orbital hyperfine fields. The latter can be 
understood as follows (see also Tab. HTjl : in orbitals with 
opposite TO, the electrons move in opposite directions. 



TABLE II: Overview whether fiorb, Bort, ^dip and Vzz depend 
on the shape of the 4f orbital (given by the absolute value of 
m), the direction of motion of the electron in that orbital 
(given by the sign of m), and whether they depend on the 
charge or on the spin of the 4f electron in a given m-orbital 





f^orb 


Borfa 




v.. 


shape of orbital 


yes 


yes 


yes 


yes 


direction of motion 


yes 


yes 


no 


no 


charge or spin 


charge 


charge 


spin 


charge 



because opposite m (z-component of the orbital angular 
momentum) mean that the angular momenta of those 
orbitals have different orientations. Hence, the orbital 
fields will be opposite as well. As a function of Z, the 
contribution due to each m-orbital increases. A linear fit 
is possible (full line). In order to verify whether this is 
accidental or not, we did the same calculations for free 
lanthanide atoms and free lanthanide 3^ ions. For the 
free ions, almost the same perfect linear correlation was 
found as for the solid (dotted line in Fig. EJ). For free 
neutral atoms, the fields were slightly larger (at most 
10% for the orbital field, and 5% for the dipolar field). 

Additional to this large 4f contribution, there is also 
a 5p contribution to the orbital HFF (Fig. (Sj-c-e). The 
up and down contributions arc quite large (40-90 T), 
but they cancel each other, yielding a negligible (< 
3 T) total contribution for the 5p-orbital HFF. This 5p- 
contribution to the orbital HFF does not depend on the 
f-configuration. Taking Dy as an example - 2 electrons in 
the unfilled spin channel - the total 5p orbital field will 
be -2 T (72 T for 5p-up, -74 T for 5p-down), irrespec- 
tive whether these 2 electrons are e.g. in the m=+3 and 
m=+2 or in the m=-l and m=0 orbitals. This will be 
different for dipolar hyperfine fields and for the electric- 
field gradient. Looking separately to the 5p up and down 
contributions, one can see that they don't vanish for La 
(4f empty) and Lu (4f full), and that they increase with Z. 
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La Ce Pr NdPmSm Eu Gd Tb Dy Ho Er TmYb Lu 



FIG. 5: Contribution of each m-orbital to the (a) 4f-up, (b) 
4f-down orbital moment and the contribution of the (c) 5p- 
up, (d) 5p-down, (e) 5p up-l-down electrons to the orbital 
magnetic moment of a lanthanide in Fe. Data points: results 
from calculations for lanthanides in Fe. Full lines: linear fit 
through these data points. Dotted lines: linear fit through a 
complete set of calculations for free lanthanide ions. 




La Ce Pr NdPmSm Eu Gd Tb Dy Ho Er TmYb Lu 



FIG. 6: Contribution of each m-orbital to the (a) 4f-up, (b) 
4f-down orbital HFF and the contribution of the (c) 5p-up, 
(d) 5p-down, (e) 5p up+down electrons to the orbital HFF 
of a lanthanide in Fe. Data points: results from calculations 
for lanthanides in Fe. Full lines: linear fit through these data 
points. Dotted lines: linear fit through a complete set of 
calculations for free lanthanide ions. 
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La Ce Pr NdPmSm Eu Gd Tb Dy Ho Er Tm Yb Lu 



FIG. 7: Contribution of each m-orbital to the (a) 4f-up, (b) 
4f-down dipolar HFF and the induced (c) 5p-up, (d) 5p-down, 
(e) 5p up-|-down contributions to the dipolar HFF of a lan- 
thanide in Fe. Data points: results from calculations for lan- 
thanides in Fe. Full lines: linear fit through these data points. 
Dotted lines: linear fit through a complete set of calculations 
for free lanthanide ions. 
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La Ce Pr NdPmSmEu Gd Tb Dy Ho ErTmYbLu 



FIG. 8: Contribution of each m-orbital to the (a) 4f-up, (b) 
4f-down to the total Yzz and the induced (c) 5p-up, (d) 5p- 
down, (e) 5p up-|-down to the total Vzz of a lanthanide in 
Fe. Data points: results from calculations for lanthanides in 
Fe. Full lines: linear fit through these data points. Dotted 
lines: linear fit through a complete set of calculations for free 
lanthanide ions. 



9 



Moreover, we observed that this contribution disappears 
if the SO-coupUng is switched off. Therefore we can con- 
clude that these 5p-contributions are due to an intrinsic 
p-effect, induced by the SO-couphng on the p-electrons 
which breaks the cubic symmetry. FinaUy there is also a 
contribution from the valence 6p-electrons, also induced 
by the SO-coupling, but this contribution is really small 
and can be neglected. The numerical stability of this 
analysis can be checked by repeating it for the orbital 
moment (Fig. The contribution to the orbital mo- 
ment of each m-orbital should be m by definition, which 
is found indeed (Fig.[31-a,b). The 5p-contributions to the 
orbital moment have the same origin as in the case of the 
orbital HFF, namely the SO-coupling. 

For the dipolar HFF as well, the 4f-contribution re- 
mains the dominant one (Fig. El-a-jb), but it is one order 
of magnitude smaller than the orbital HFF. The system- 
atics are different from the orbital moment and orbital 
HFF as well (Fig. |7Ja,b and Tab. HTJ. First of aU, the 
dipolar HFF does not depend on the direction of motion 
of an electron, such that ±m-orbitals yield the same dipo- 
lar HFF. Secondly, Bdip depends explicitly on the elec- 
tron spin, such that an electron with opposite spin in the 
same m-orbital yields an oppositie field. Furthermore, we 
can observe from Fig.El-c,d,e that the 5p-contributions to 
Bdip depend on the 4f-occupation. If we take again Dy 
as an example and put the 2 electrons in the -1-3 and +2 
orbitals (spin up) we get an induced 5p-contribution of 
4 T (-8 T for 5p-up and 12 T for 5p-down), while this 
is -3 T (4 T for 5p-up and -7 T for 5p-down) if the 2 
electrons are in the +1 and -2 orbitals (spin up). Such 
a 4f-dependence was not present for fiorb and Borb- In 
Sec. lIVBl we will see that also for the EFG there is such 
an explicit 4f-dependence, and we will be able to explain 
this by the radial dependencies, which are 1/r for iiorb 
and Borb , and 1 /r^ for Bdip and Vzz ■ Another obser- 
vation from Fig. is that the induced 5p-contributions 
behave differently in the first and the second half of the 
lanthanide series. This suggest a spin-dependent inter- 
action: In the first half of the series the unfilled f-band 
is the down band. The 5p-up contribution is large, the 
5p-down is smaller. In the second half the unfilled f-band 
is the up band. Now the 5p-up contribution is small and 
5p-down smaller. Apparently the 4f-electrons induce a 
larger contribution in the 5p with opposite spin. The 6p- 
contribution remains negligible also for the dipolar HFF. 

For the Fermi contribution, we cannot follow the pro- 
cedure of Fig. El as this contribution is less directly con- 
nected to the occupation of the f-states. In Fig. O the 
Fermi contribution (for valence (mainly 5s and 6s) and 
core electrons (mainly Is to 4s) separately) is plotted as 
a function of the 4f spin moment, and this for different 
lanthanides and different types of solutions. Both the 
core and valence Fermi contribution depend linearly on 
the 4f moment, and they sum to an almost constant value 
that is independent of the 4f moment. As this constant 
contribution is there even for a zero 4f spin moment, it 
must be due to other than the 4f electrons. Actually, it 
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FIG. 9: Fermi contact contribution to the HFF as a function 
of 4f spin moment for a large set of different solutions for 
different lanthanides in Fe. White circles: core contribution. 
Black circles: valence contribution. Gray circles: total Fermi 
contact contribution. The lines through the core and valence 
contributions are linear fits, the line through the total Fermi 
field is the sum of those two fits. 



results from the 6s-polarization due to the small and rel- 
atively constant 5d-moment. The 4f spin moment is spa- 
tially outside the core s-orbitals, but inside the valence 
s-orbitals. It will therefore induce exactly the opposite 
s-polarization on both of themS^. 

Encouraged by the above observation that the m- 
orbitals of lanthanides in Fe yield identical orbital and 
dipolar fields as for free lanthanide ions, we now suppose 
that the ground state occupation of the 4f shell for lan- 
thanides in Fe is the same Hund's rules ground state as 
for free lanthanide 3"'' ions. 

For Tm"^+ as an example, Hund's rules prescribe the 
orbitals with m=+3 and m=+2 for the 4f-up spin to 
be unoccupied, which yields from Fig. and Fig. [7| the 
orbital and dipolar contributions listed in Tab. IIIII The 
Fermi contact field is almost independent on the orbital 
occupation, and is -39 T for Tm and all other lanthanides 
in Fe. The total HFF for the ground state of Tm in Fe - 
corresponding to a Tm^+ configuration - is hence -745- 
53-39=-837 T. This value is very close to what was found 
in a direct way in 'case 1' of Tab.Ql a solution that has 
the same type of occupation as the Hund's rules ground 
state. 

In this way, we can obtain by Fig. an LDA-I-U or- 
bital HFF for all lanthanides - in the antiferromagnetic 
type of solution - which is given by the full line in Fig.|2|- 
a. In the same way. Fig. [T] leads to the LDA-I-U dipolar 
field, given by the full line in Fig. and Fig. to the 
LDA-fU Fermi field in Fig.Ql-c. (It should be noted here 
that the core part from this calculated Fermi field most 
likely suffers from the typical 'LDA core error', for which 
recently a promising cure has been proposedS&.) The sum 
of all these contributions leads to an LDA-I-U value for 
the total HFF of all lanthanide 3+ ions in Fe, which is 
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TABLE III: Contributions to Borb, Bdip and Vzz for Tm in 
Fe in the Hund's rules ground state, using the information 
from figs. I6I8I 



Tm 


Bor6(T) 


Bdip(T) 


V,,(102iV/m2) 


4f-up 


-745 


-49 


-76.2 


4f-down 


-1 


-1 


0.9 


5p-up 


-84 


10 


16.1 


5p-down 


83 


-13 


21.2 


6p-up 


10 





-0.7 


6p-down 


-8 





0.1 


Sum 


-745 


-53 


-38.6 



given by the full line in Fig. ^b. The agreement with 
experiment is rather good for many lanthanides, espe- 
cially in the middle of the series. Towards the beginning 
and the end, there is moderate disagreement. Only for 
Yb there is a large deviation between an accurate and 
reliable experimental value of -125 T (TDPAC) and a 
calculated value for the ion of -570 T. Yb, however, 
often occurs in a divalent configuration, where the 4f shell 
is completely filled. This results in the -39 T of the Fermi 
contact field only, which is in much better agreement with 
experiment. Nevertheless, a trivalent configuration was 
suggested befor from the following experimental con- 
siderations: the -125 T for Yb in Fe was considered to 
be 'large', much larger than the (-)61 T for the divalent 
Lu in Fe. This additional -64 T was taken as stemming 
from an orbital contribution, which must lead to the con- 
clusion that the Yb is trivalent. From Fig.^a, however, 
we see that such a difference of 64 T is almost negligible, 
and that the expected orbital contribution for a triva- 
lent state would be 10 times larger. The HFF for the 
other lanthanides in a divalent configuration is shown in 
Fig. ^b as well (dashed line). Although the agreement 
is better now for Er and Tm, the uncertainty on the cal- 
culated values is too large to pretend that they would be 
divalent as well (which is for Er extremely unlikely, any- 
way). Another lanthanidc that is often divalent is Eu. 
The experimental HFF lies halfway between the HFF for 
Eu^+ and Eu"^"*", such that at this stage no definite con- 
clusion on the valency is possible yet. Mossbauer isomer 
shift data™^, however, convincingly point to trivalency 
for Eu in Fe (see Sec. lIVBl for a continuation of this dis- 
cussion and a final conclusion on the Eu valency). Irre- 
spective of the valence configuration we consider, one can 
see that a moderate deviation from experiment for La, Ce 
and Pr still remains for the LDA+U calculations. The 
LDA results are closer to experiment. As was mentioned 
in Sec. LDA calculations for lanthanides represent an 
itinerant (also called delocalized) 4f configuration, which 
is mostly not what is found in Nature: the radius of the 
4f orbitals is not very large, overlap with the orbitals of 
neighboring atoms is negligible, and as a result the 4f 
orbitals are localized. The only exception is Ce. As the 



4f radius gets smaller for increasing atomic number Z, 
the 4f orbitals of Ce reach most outwards. The radius 
is just large enough to allow overlap and hence delocal- 
ization in materials where the nearest neighbor distance 
is not too large. Therefore, Ce can be either trivalent 
or itinerant, depending on the materia&Si. One can 
say that in the lanthanidc series there is a delocalization- 
localization transition, that happens already at the very 
first element, Ce. The good agreement between experi- 
ments and the LDA-results up to Pr in Fig.QJb, suggests 
that not only Ce but also Pr has delocalized 4f electrons 
in an Fe host: a 'postponed' delocalization-localization 
transition. Three additional arguments support this hy- 
pothesis. The first one is given in Fig. ^1 where for all 
lanthanides in Fe the number of electrons inside the muf- 
fin tin sphere with radius Rmt = 2.45 a.u. is plot for the 
3 most occupied m-orbitals in the unfilled spin channel, 
and compared with the average number of electrons in 
the same muffin tin sphere for occupied orbitals in the 
free lanthanidc ions (which is almost unity). All data are 
for LDA-|~U calculations. Even though LDA-I-U favors 
localization, one can see that the occupation of the or- 
bitals for Ce and Pr is far from 1: the 4f electrons are 
distributed in noninteger quantities over all 7 m-orbitals 
and/or reach out of the muffin tin sphere. This is a sig- 
nature of delocalized behavior. (Two notes: 1. our pro- 
cedure to extract the contributions for a single orbital 
works even with such noninteger occupations as input. 
The LDA-f-U values in Figs, ^b and ^] are therefore 
really values for localized orbitals, although this local- 
ization is not realized in the LDA-I-U calculations them- 
selves; 2. Fig. ^] should not be considered as a proof 
that from Nd on the 4f electrons are localized. These are 
results from LDA-I-U, which favors localization. Even 
LDA-I-U is not able to localize the 4f up to Pr. Therefore 
fig. ^] offer a lower bound for the delocalization transi- 
tion (delocalization at least up to Pr). It does not ex- 
clude that, in Nature, Nd, Pm,... in Fe could be delocal- 
ized as well.) A second supporting argument comes from 
the lanthanide-Fe distance: 2.60 A, which is considerably 
smaller than a typical lanthanide-lanthanide distance in 
pure lanthanidc metals (4.08 A). One can expect from 
this a large overlap between the 4f orbitals and the Fe- 
3d, and hence a stronger tendency to delocalization. This 
argument will be further quantified in terms of pressure 
in Sec. lV AI Finally, in Sec. lIVBl we will show that also 
the EFG of Ce in Fe indicates delocalized 4f behavior. In 
conclusion, we can say that there are several strong in- 
dications that the delocalization-localization transitions 
for lanthanides in Fe is somewhat postponed, at least up 
to Pr. We will come back to this in Sec. IV AI 

Regardless of the ambiguity to choose between diva- 
lent and trivalent configurations, one conclusion can be 
made unambiguously from Fig. ^b: the magnitudes of 
the HFF's with LDA-I-U are much closer to experiment 
than with LDA, certainly for the heavier lanthanides. 
Together with the antiferromagnetic coupling which is 
reproduced by LDA-I-U but not everywhere by LDA 
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troscopy from Ref. |43. Niesen and Ofer. however, have 
later shown^ by ^^^Eu Mossbauer spectroscopy that 
14^ = -11.110^1 V/m^. For Ce m Fe the quadrupole 
coupHng constant (which contains the product between 
Vzz and the quadrupole moment Q) has been determined 
by ^"'^Ce NMRS^ to be almost zero. The quadrupole 
moment for ^'^^Ce is not known, but assuming a typi- 
cal value of 1 barn leads to a V^z that is practically zero 
(gray symbol in Fig. Ill() . The zero Vzz for Yb is not ex- 
plicitly mentioned in the literature, but can be inferred 
from Fig. 1 in Ref. 112, which shows a purely magnetic 
interaction. 



La Ce Pr NdPmSmEu Gd Tb Dy Ho Er Tm Yb Lu 



FIG. 10: The average occupation of the m-orbitals for the 
free ions (full line) compared with the occupation of the 3 
orbitals in the unfilled spin-channel that are most occupied 
for lanthanide in Fe (dotted, dashed and dotted-dashed lines) 
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FIG. 11: Vzz for lanthanides in Fe. Symbols: experimental 
datai-i— (see text). The gray symbol for Ce indicates 
that this value is a - rather safe - guess (see text). Dotted- 
dashed line: LDA results. Full line: LDA-I-U values for triva- 
lent lanthanides (using the information of Fig. |S1 see text). 
Dashed line: LDA-fU values for divalent lanthanides. 



(Fig. 21) . this is hard evidence for the fact that LDA-I-U 
performs considerably better than LDA also in these sys- 
tems. 



IV. ELECTRIC-FIELD GRADIENTS 

A. Experimental data set 

Only few experimental data on the main component 
Vzz of the electric-field gradient tensor for lanthanides 
in Fe are available (Fig. Ill|l . For Eu, Dy, Er and 
Tm, Mossbauer measurementsiS*^ were done. Some 
attention is needed for Eu, which is reported to have 
Vzz = in Ref. IT^ based on ^^^Eu Mossbauer spec- 



B. LDA+U calculations 

As in an LDA-calculation all f-orbitals are roughly 
equally populated, there is almost no spatial anisotropy 
and Vzz will be close to zero (dotted-dashed line in 
Fig. Ill|l . We therefore turn immediately to LDA-I-U cal- 
culations. We are confronted with the same problems as 
described in Sec. IIII CI - Vzz sensitively depends on the 
f-elcctron density matrix and there is no criterion to de- 
termine which density matrix corresponds to the ground 
state - and therefore we determine again the individual 
contribution to Vzz of every 4f m-orbital and the induced 
5p contributions to Vzz (Fig. (Sjand Tab. ITT|) . Just as for 
the dipolar HFF, the EFG due to the 4f orbitals them- 
selves does not depend on the direction of motion of the 
electron, and ztm-orbitals yield the same V^^z- On the 
other hand; the EFG depends on the charge and not on 
the spin, such that also up and down 4f-electrons yield 
the same Vzz- The 5p contribution depends on the 4f 
occupation, just as for the dipolar HFF. Assuming a 
Hund's rules type of occupation either for divalent or 
trivalent ions, leads to reasonable agreement with exper- 
iment fFig. lll|) . Fig. ^2 provides further evidence for the 
fact that Yb in Fe really is divalent: trivalent Yb has 
a very large Vzz, while the experimental value is zero, in 
agreement with the divalent prediction. For Eu, the mea- 
sured negative ]4z is in good agreement with the calcu- 
lated value for a trivalent state. Together with the exper- 
imental evidence based on the isomer shift fSec. Illimi . 
we can now firmly conclude that Eu in Fe is indeed triva- 
lent, and has a valence state which is different from the 
one of Yb. For Er, Fig. 1111 suggests trivalency as well, a 
conclusion that could not be unambiguously made based 
on the HFF only fSec. IIII^ . For Tm, the situation re- 
mains undecided. Irrespectively of the valence state we 
consider for Ce, a large deviation from experiment can 
be observed. However, the experiment matches very well 
with the itinerant LDA result. This is another indication 
that Ce in Fe has delocalized f-electron. 

Now we can analyze which are the main electrons that 
provide the anisotropy that leads to the EFG. It has been 
shown beforeS^ in a rigorous way that for metals with 
S-, p- and d-electrons the total Vzz can be obtained as 
a sum of a quantity V^z~^ Kz""^ (neglecting small 
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TABLE IV: Contributions to Kz for Tb in Fe, with the down- 
channel for Tb-4f completely filled and with for the 4f-up 
channel one electron in the m=-l orbital. The first column 
gives the rigorous notation of each contribution (see R.ef . l6fl^ . 
In the second column an interpretative notation is defined, 
that is used in the text and in Fig. |H| The third column 
gives the energy region in the Density Of States (DOS) near 
to which these states are found (Ep means "near the Fermi 
energy", negative values are below the Fermi energy). Units; 
10^1 V/m2. 







in DOS (eV) 


up 


down 


V^-^ (5d) 




Ep 


0.1 


0.4 


Vi-^ (5p) 




-23 


5.6 


9.6 


V^-^ (6p) 


^ zz 


Ep 


-0.2 


0.4 


V!-' (4f) 


^4f 
^ zz 


-5 


-29.0 


2.5 


14z 






-23.5 


12.9 sum = -10.6 



contributions from the interstitial region of the crystal). 
They measure the nonspherical p and d charge densities 
(^^{t) and p^^{r\ respectively, weighted by an inte- 
gral over 1 jr^ : 



rR „P-P 



-P cc f 
Jo 



.R d-d 



(1) 

(2) 



R is the radius of the muffin tin sphere of the considered 
atom. The factor strongly emphasizes the contribu- 
tion from the region close to the nucleus, with small r. 
s-clcctrons do not contribute as they have spherical sym- 
metry, and so-called 'mixed' s-d or s-p contributions are 
negligible and therefor omitted. This can be extended to 
materials with f-electrons, such that V^z for lanthanides 
can be written as: 



Vz, 



yP-P + yd-d , yi- 

* 7.7. ' * 7.7. ' * 7.7. 



(3) 



We now apply this analysis to Tb in Fc, which is a par- 
ticularly clear example because all 4f-down orbitals are 
fully occupied and there is only a single 4f-up electron. 
Tab. IIVI shows the different contributions to Vz^ when 
this single electron is put in the m=-l (up) orbital (this 
is not the ground state, but this is just an example, any- 
way). Tab. II Vl shows that the main contribution to the 
total Vzz = -10.610^1 Y/ia is due to the single 4f-up 
electron: V^^J = -29.0. There is a large contribution of 
5.6+9.6=15.2 with the opposite sign due to p-electrons. 
What is surprising is that this p-contribution does not 
stem from the valence 6p electrons, but from the entirely 
filled and strongly bound 5p shell, which lies more than 
20 eV below the Fermi energy. Intuitively, one would 
have assumed such a filled and well-bound shell to be en- 
tirely spherically symmetric, which would mean V^^ ~ 0. 
And indeed, the 5p-anisotropy Ap = i (rip^ -I- ?T-py) ^np^ 
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FIG. 12: a) The anisotropic p-density p2o'^{''') for Tb in Fe 
(sum of up and down electrons, arbitrary units for y-scale). 
b) Left axis (arbitrary units): p2o^{'r)/r^ for Tb in Fe. Right 
axis (arbitrary units): integral of p2o^{''')/r^, which is called 
Kz ^^ (apart from a constant factor with negative sign), c) 
and d): idem, but for the f-f contribution. 



is very small: 0.0030 (up) and 0.0040 (down) (it is shown 
in Ref. (i!) that Ap is proportional to V^J' ; Up^ is the 
number of electrons in the pi orbital, and Ap measures 
the unequal occupation of the 3 p-orbitals). However, a 
considerable part of this anisotropy stems from a region 
very close to the nucleus and hence gets amplified by 
the factor. This is demonstrated in Fig. ^1 where 
for the same Tb-configuration the bare anisotropic p-p 
charge density P20~'^(^) shown, before (Fig. 1121-a) and 
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after (Fig.^J-b) weighing with a 1/r^ factor, and also af- 
ter integration fFig.lT^b. right axis). Figs.lT^c andlT^d 
repeat this for the f-f contribution. The final integrals are 
clearly determined exclusively by anisotropics in a region 
closer than 0.05 A to the nucleus. In the expression for 
a spin dipolar field, the same factor Ijr^ is present, ex- 
plaining why a similar dependence of the 5p dipolar field 
on the 4f occupation was observed there (Fig. [71-c,d;e). 

A generalization of this analysis is given in Fig. |H1 
which shows the same individual 4f m-orbital contribu- 
tions to for lanthanides split into 4f and 5p (i.e. when 
a given 4f m-orbital is occupied, Fig. |Hl-a,b shows the di- 
rect contribution from this orbital, while Fig. |Hl-e shows 
the corresponding induced contribution of the 5p-shell 
(up and down summed)). This shows that the oppo- 
site signs for 5p and 4f as seen in the example of Tb is 
a general effect: occupying the 4f m=0 orbital gives a 
negative direct 4f contribution but induces a positive 5p 
contribution, etc. This can be understood as follows. A 
negative T4z corresponds to charge accumulation along 
the z-axis, a positive 14z to charge accumulation in the 
xy-plane. The shape of the 4f orbitals is such that m=0 
has its charge mainly along the z-axis (reflected in a neg- 
ative 14t^), while the xy-plane is more and more occupied 
for larger \m\. Apparently the 4f electrons dispel the 5p 
electrons: if m=0 is occupied, then the 5p electrons are 
forced away from the z-axis into the xy-plane, resulting 
in a positive V^^ and Ap. As Tab. IIVI shows, a 4f-up 
electron distorts the 5p orbitals with either spin: this is 
an interaction between the electron charges, not between 
the spins. 

Recently, the EFG of the actinide U has been analyzed 
in UO2 by R. Laskowski et al. fRef . ItoI) . using the same 
APW-l-lo method as used in this work. These authors 
show in their Fig. 2 the contributions of V^^'^, Kz~^ ''^^'^ 
V^~^ as a function of the deformation of the oxygen cage 
that surrounds the U atom. They do not further divide 
the p-p contribution in 6p and 7p (here the 6p shell is 
entirely filled and well-bound). Without deformation of 
the oxygen cage, the crystallographic surrounding is cu- 
bic and the slightly non-zero 14z ~ —2 10^^ V/m'^ is due 
to spin-orbit coupling only, as is the case for lanthanides 
in Fe (compare with Nd in Fig. Ill|l . This small 14z is a 
sum of a f-f contribution of -1-20 and a p-p contribution 
of -21 (the d-d contribution is small: -1). This can be 
compared to Fig. |S1 for Nd with the m=(-3, -2, -1) or- 
bitals filled, which leads to a 4f-contribution +14 and a 
5p-contribution of -12, quite similar values. Because of 
this analogy, we suggest that also for U in UO2 th p- 
p contribution for an undistorted oxygen cage is due to 
the completely filled 6p shell. This interpretation would 
furthermore imply that as a function of oxygen cage de- 
formation, the 6p-contribution in Fig. 2 of Ref. [t^ would 
remain almost constant (just as the 5f-contribution does), 
and that the strong decrease of the total p-p contribution 
is due to 7p only. This makes sense, as this decrease is 
attributecffi to the tails of the 0-2p wave functions, and 
hence should appear near the Fermi energy (= the region 
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FIG. 13: The experimental transition pressure from localized 
4f to delocalized 4f electrons for pure lanthanides (squares) 
and two possible choices for the effective pressure felt by lan- 
thanides in Fe (full and dotted lines) (for details see text) 



of the 7p). 



ELABORATIONS 



Pressure 



It is well known experimentally^, that the pure lan- 
thanides exhibit a large variety of structural phase tran- 
sitions as a function of external pressure. At certain pres- 
sures, volume collapses are sometimes observed and at- 
tributed to the delocalization of the f electrons. These 
delocalization pressures have been determined for 6 el- 
ements of the lanthanide series: CeSSiZ^, pj^74,75,76,77^ 
N^jT^^g^^ Srrmi^, Gd2i and Dy^ (Fig.^. For Ce 
the delocalization of the f electrons occurs around the 
pressure of 1 GPa and is accompanied by a volume col- 
lapse of 16% at the isostructural transition to another 
fee phase (a-Ce). Pr transforms to a a-U structure at 
20 GPa with a volume collapse of 9-12%. In Gd 4f delo- 
calization occurs at 59 GPa when the structure changes 
to body-centered monoclinic (bcm) with a volume col- 
lapse of 10.6% and in Dy this happens at 73 GPa with 
a volume collapse of 6%. For Nd and Sm no volume 
collapse has been observed. In these two cases the delo- 
calization of the 4f electrons was associated with the ap- 
pearance of low-symmetry structures (similar with those 
that appear in Pr, Gd and Dy cases with volume col- 
lapse). This is a somewhat ambiguous procedure. For 
Nd two transition pressures have been proposed: 40 GPa 
(corresponding to the transition to an hP3 structuraSiS) 
and 113 GPa (corresponding to the transition to the a-U 
phased) . For Sm the delocalization pressure is proposed 
to be 37 GPa, when the Sm structure changes to hP3. 
For Pm we have only a lower limit for the transition 
pressure, 60 GPa (until this pressure no low-symmetry 
structure has been observed^). In our calculations for 
lanthanide impurities we replace an Fe atom from an iron 
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lattice by a lanthanidc atom. Obviously the lanthanide 
atom (which has a much larger volume) will feel a chem- 
ical pressure or effective pressure. How large will this 
effective pressure be? We concluded in Sec. IIII CI that 
at least Ce and Pr are delocalized. Hence, the effective 
pressure - which we assume to be independent on the 
lanthanidc in a first approximation - should be at least 
20 GPa (Fig. 1131). The hypcrfinc fields for Nd and Sm 
are only very approximately measured (Fig. ^a), such 
that one cannot conclude whether they are localized or 
not. Based on isomer shift and EFG (Fig. [TT)l . Eu is 
definitely localized, as are all heavier lanthanides. There- 
fore, two qualitatively different proposals for the effective 
pressure are possible: about 25 GPa (everything starting 
with Nd is localized) or about 40 GPa (everything be- 
low Eu is delocalized, except for Pm and maybe Nd). 
Assuming an effective pressure that is not constant (mo- 
tivated by the decreasing volume of heavier lanthanides) 
does not change this picture (dotted lines in Fig. ^1 - 
these lines qualitatively take the lanthanide contraction 
into account). A more accurate experimental determina- 
tion of HFF and EFG for Nd, Pm and Sm in Fe would 
allow to distinguish between both scenarios and would 
allow to determine the real position of the delocalization- 
localization transition in this system. Experimental data 
for Nd in Fe are also for another reason interesting: if 
Nd in Fe would be found to be localized (itinerant) and 
the effective pressure of 40 GPa would known to be cor- 
rect (from a Sm-measuremcnt for instance) than the dc- 
localization pressure of 113 GPa (35 GPa) for bulk Nd 
is probably correct. If the effective pressure of 25 GPa 
would be correct, no such conclusion can be made. 

In conclusion for this section, accurate measurement 
of HFF and EFG for Nd, Pm or Sm in Fc would offer a 
lot of information. 



B. Free lanthanide ions 

LDA+U calculations for free lanthanide 3"*^ ions were 
already mentioned in Figs.|Sl-|Hl A quite complete exper- 
imental data set exists for this situation as well, both for 
the HFF (Fig. Ula) and the EFG (Fig. Ulb) (data are 
copied from Ref. l85l the original data are in Ref . l86|) . In 
the experiments, the ions were not really free but were 
incorporated in a paramagnetic salt, and the effect of 
crystal fields was removed later in order to find the free 
ion values. By the same procedure as explained for lan- 
thanides in Fe, we extracted LDA+U predictions for the 
HFF and Vzz for divalent and trivalent lanthanide ions, 
which arc given by dashed and full lines in Fig. 1141 re- 
spectively. For lanthanides in Fe, the positive z-direction 
was naturally defined by the moments of the ferromag- 
netic Fe host atoms. For free lanthanides, the total 
(=spin-|-orbital) angular momentum J determines the 
positive z-direction. Due to this different choice of axes, 
there is an apparent sign change for the heavy lanthanides 
between Fig. ^ (in Fc) and Fig. 1141 -a (free) . The agree- 
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FIG. 14: a) Experimental value for the HFF in free lanthanide 
ions, compared with LDA-I-U predictions based on Figs.|Sl|7| 
and|^ both for divalent and trivalent lanthanides. b) Exper- 
imental value for Kz in free lanthanide ions, compared with 
LDA+U predictions based on Fig. |H| both for divalent and 
trivalent lanthanides. 



mcnt with experiment is again quite nice. Eu is divalent 
(it was trivalent in Fc), while Yb is trivalent (it was di- 
valent in Fe) . For Sm there is a large deviation, both for 
the HFF and V^z. But this is no surprise: it is well-known 
there are low-lying excited states in Sm which will mix 
with the ground state, such that our procedure which is 
based on Hund's rules ground states is expected to fail. 



C. Non-collinear Magnetism 

The possibility to consider non-collinear magnetism 
at every infinitesimal region of space has recently been 
implemented^^ in the WIEN2k code. In principal this 
can be an important feature even for collinear antiferro- 
magncts as we are dealing with here: it allows the spin 
moment to turn gradually from the Fe-oricntation to the 
opposite lanthanide-orientation, and this is a better repli- 
cation of what happens also in nature. We did not at- 
tempt a full study, but calculated HFF and EFG for Tm 
in Fe only. All technical parameters were chosen exactly 
the same as in the collinear calculations. For Tm in Fe as 
test example, the total HFF changes from -822 T in the 
collinear LDA+U calculation to -838 T in a non-collinear 
one, while the EFG remains exactly the same: -38.1 lO'^^ 
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V/m^. Such a change of 16 T is not small in absolute 
value, but is negligible compared to the large values of 
the HFF's in this problem. Therefore we conclude that 
non-coUinear magnetism doesn't play an important role 
for lanthanides in Fe. 



VI. CONCLUSIONS AND OUTLOOK 

We have demonstrated that LDA leads to a qualita- 
tively wrong behavior for the magnetism of lanthanides 
in Fe. Using LDA+U, qualitative and quite reason- 
able quantitative agreement with experiment is obtained, 
both for HFF and EFG. This shows that the semi ab ini- 
tio LDA+U method is a useful tool even for such sensitive 
quantities as hyperfine parameters of strongly correlated 
impurities in an itinerant magnetic host. We could come 
to these conclusions only after applying a careful strat- 
egy in order to cope with the lack of a good criterion 
to determine the true ground state if LDA+U is used. 
For all lanthanides the 4f spin moment couples antifer- 
romagnetically to the Fe 3d moment, in agreement with 
the model of Campbell and Brooks. The orbital HFF 
is by far the dominant contribution to the total HFF 
(Figs, ^andinj. The calculated value of the EFG for 
lanthanides in Fe agrees well with the few experimental 
data (Fig. [TT|l . We discovered an unexpectedly strong 
contribution of the completely filled 5p shell to the dipo- 
lar HFF and the EFG, which can be explained by their 
common dependence: small deformations of the 5p 
shell in a region close to the nucleus are strongly empha- 
sized. A rcinterprctation of recent EFG calculationsSi 
for Uranium in UO2 suggests that the same is true for 
the 6p shell in actinides. Furthermore, we conclude that 
Yb is divalent in an Fe host, while all other lanthanides 
are trivalent (including Eu, with perhaps an exception 
for Tm). The lightest lanthanides (at least up to Pr) 
show delocalized 4f behavior, and we conclude that the 
delocalization-localization transition that typically hap- 
pens already at Ce is postponed for lanthanides in Fe: it 
falls at least after Pr and current experiments do not ex- 
clude that it could go up to Sm (although Pm is certainly 
localized). This can be explained by the large effective 
pressure that is felt by these lanthanidc impurities (ei- 
ther 25 or 40 GPa - Fig. I13|l . leading to a larger overlap 
between the 4f wave functions and the neighboring Fe- 
3d. The question of a postponed localization transition 
has never been touched before in the 40 years of exper- 
iments on this system. This illustrates what can be the 
added value of ab-initio calculations for hyperfine inter- 
actions studies. Also in the case of free lanthanide ions, 
HFF and EFG can be quantitatively reproduced. Re- 
markably, Eu is divalent in this case and Yb is trivalent 



- just the opposite as for lanthanides in Fe. The effect of 
fully non-coUinear magnetism on this problem was tested 
to be negligible. 

Obviously, the numerical agreement of the calculated 
hyperfine fields with experiment is for these lanthanides 
still much less satisfactorily then it is for lighter impuri- 
ties. With the current methods, there arc only limited 
possibilities to improve the accuracy of the calculations. 
The supercell can be extended to e.g. 32 atoms and re- 
laxation of the Fe-neighbours can be calculated for every 
individual element (this requires the calculation of forces 
including spin-orbit coupling and LDA+U, which is time- 
consuming and not yet fully implemented in WIEN2k). 
But as inevitably a rather arbitrary choice remains to be 
made for the value of U, it is not clear whether these 
sophistications will really improve the agreement with 
experiment. And most likely they will not add anything 
new to the physical insight. In our opinion, new progress 
in this topic will have to come from experiment. Many of 
the experimentally determined HFF's and EFG's carry 
still large error bars. Accurate measurements - for in- 
stance with the NMR/ON method - are desirable (note 
that NMR/ON has not yet been applied for any of the 
lanthanides with a large EFG and/or HFF: such large 
hyperfine interactions put severe requirements on the 
equipment). The predicted HFF's and EFG's from this 
work should allow to reduce considerably the frequency 
domain that has to be scanned in an NMR/ON experi- 
ment, and warrants a more physical and reliable interpre- 
tation of the observed resonances. As most worthwhile 
experiments, we suggest a more accurate determination 
of HFF and EFG for Pr to Sm: this would allow to ex- 
amine experimentally the position of the delocalization- 
localization transition. Once a data set with improved 
accuracy will be available, it can serve in its turn as a 
testing ground for future generations of ab initio many 
body methods. 
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